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Foreword 


This  report  was  prepared  by  Urho  A.  Uotila  and  Richard  H.  Rapp, 
Professors,  Department  of  Geodetic  Science  at  The  Ohio  State  University, 
under  Air  Force  Contract  No.  F19628-76-C-0Cn0,  OSURF  Projects  No. 
710334  and  710335.  This  is  a final  report  of  the  contract  covering  time 
period  July  1,  1975  to  September  30,  19  78.  It  has  been  administered  by 
the  Air  Force  Geophysics  Laboratory,  Air  Force  Systems  Command, 
Hanscom  AFB,  Massachusetts  with  Mr.  Bela  Szabo,  Contract  Monitor. 

The  research  done  under  this  contract  has  been  previously  reported 
in  twenty-two  scientific  and  two  internal  reports.  In  the  following  the 
scientific  work  done  under  this  contract  will  be  summarized.  The  sub- 
headings 1.1-1.  7,  and  2,  are  written  by  Uotila  and  1.8-1.12,  3,  4,  and  5 
are  written  by  Rapp  related  to  research  work  done  under  OSURF  projects 
710334  and  710335,  respectively,  which  have  been  under  their  supervision. 

The  authors  wish  to  thank  all  of  those  who  have  participated  in  the 
research  under  the  contract.  Special  mention  should  be  made  about  the 
excellent  contributions  of  Drs.  Hajela,  Jekeli,  Kearsley,  Moritz,  Rapp, 
Rummel,  Schwarz,  Sjoberg,  Sunkel  and  Tscherning  to  the  successful 
completion  of  the  research  contract.  The  authors  acknowledge  the  cooper- 
ation and  support  given  and  express  their  appreciation  to  the  Contract 
Monitor,  Mr.  Bela  Szabo,  for  the  stimulating  technical  and  scientific 
discussions. 
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1.  Collocation  and  Related  Subjects 

Dr.  Moritz  prepared  six,  Dr.  Sjoberg  two,  and  Drs.  Hummel,  Siinkel, 
Tscherning  and  Jekeli  one  technical  report  in  this  area.  Part  of  the 
summaries  reported  under  1. 1-1.7  are  taken  directly  from  Moritz's  and 
Siinkel's  reports  without  quotation  marks.  The  titles  of  the  reports  are  used 
as  the  subheadings  for  the  presentation. 

1.1  Integral  Formulas  and  Collocation 

Various  aspects  of  the  interplay  between  least-squares  collocation  and 
classical  integral  formulas  were  subjects  of  Moritz's  report  (1975)  "Integral 
Formulas  and  Collocation.  " 

The  first  three  sections  of  the  report  were  concerned  with  the  proof  that 
integral  formulas,  such  as  Stokes’  and  Molodensky's  equations,  might  be  con- 
sidered as  limiting  eases  of  collocation  for  homogeneous  and  regularly  and 
densely  distributed  data.  Moritz  showed  that  the  two  types  of  methods  are  in- 
deed compatible  and,  what  is  more,  they  even  complement  each  other  from  the 
point  of  view  of  practical  application. 

Sections  4 and  5 dealt  with  applications  of  least-squares  methods  to  adjust 
continuous  data,  with  a view  to  using  them  in  Stokes'  and  other  integral  formulas. 

Practical  and  numerical  aspects  were  presented  in  the  last  section  and  the 
relative  merits  of  both  types  of  techniques  and  their  interplay  were  considered. 
With  respect  to  numerical  computation,  Moritz  concluded  that  integral  formulas 
and  collocation  techniques  mutually  complement  each  other,  so  that  in  many 
practical  cases  a judicious  combination  of  the  two  procedures  may  be  practically 
most  convenient. 

Throughout  his  report,  the  so-called  spherical  approximation  was  user!.  It 
consisted  in  neglecting,  in  equations  that  relate  quantities  of  the  anomalous  gravity 
field,  small  terms  on  the  order  of  the  flattening;  thereby  ellipsoidal  relations 
were  formally  transformed  into  spherical  formulas.  More  precisely,  a point 
on  or  near  the  earth's  surface,  having  geodetic  coordinates  4,  A and  height  h 
above  the  ellipsoid,  was  mapped  into  a point  that  had  spherical  coordinate's  4 , 

X and  height  h above  a sphere  of  radius  R,  R - 6370  km  being  the  mean  radius 
of  the  earth  and  4 , A , h being  numerically  the  same  in  lx)th  cases. 

The  spherical  approximation  considerably  simplified  formulas,  while  its 
error  was  negligible  in  most  practical  cases.  It  underlies  Stokes'  and  Vening 
Meinesz'  integral  formulas,  as  well  as  most  solutions  of  Molodensky's  problem; 
it  also  underlies  least-squares  collocation.  If  necessary,  ellipsoidal  effects 
could  easily  be  taken  into  account  by  small  correction  terms  to  the  spherical 
approximation  as  shown  by  Moritz  in  his  earlier  ro|)ort  (1974). 
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1.2  Covariance  Functions  in  I .cast -Squares  Collocation 


Least-squares  collocation  depends  essentially  on  the  covariance  functions 
used.  This  is  especially  true  for  accuracy  studies,  for  which,  on  the  other  hand, 
collocation  provides  a powerful  mathematical  apparatus. 

For  instance,  it  is  known  that  some  analytical  expressions  for  covariance 
functions  lead  to  imaginary  standard  errors.  What  is  wrong  with  such  functions? 
They  are  not  positive  definite. 

It  is  also  well  known  that  the  covariance  functions  are  responsible  for  the 
precise  mathematical  structure  of  the  gravity  field  through  covariance  propaga- 
tion. This  implies  that  the  basic  covariance  functions  must  be  harmonic. 

It  appears,  therefore,  appropriate  to  elaborate,  in  some  detail  mathematical 
properties  of  covariance  functions  such  as  positive  definiteness  and  harmonicity. 

Another  question  is  how  to  characterize  a covariance  function  sufficiently 
well  by  a small  number  of  parameters,  in  such  a way  that  two  different  covariance 
functions  that  have  these  parameters  in  common,  give  approximately  the  same 
result. 

To  find  a good  analytical  covariance  function,  one  tries  to  represent  it  as 
a linear  combination  of  simpler  functions.  It  is,  therefore,  desirable  to  know 
the  behavior  of  such  simple  models  which  may  serve  as  building  blocks  for  a 
global  covariance  function. 

All  these  problems  were  considered  in  the  first  part  of  Moritz's  report 
(1976a):  "Covariance  Functions  in  Least-Squares  Collocation.  " He  dealt  with 
the  mathematical  structure  of  covariance  functions.  The  properties  of  isotropy, 
harmonicity  and  positive  definiteness  were  discussed,  and  Moritz  suggested  that 
a covariance  function  may  be  characterized  by  three  essential  parameters;  the 
variance,  the  correlation  length  and  a curvature  parameter.  He  also  considered 
some  spatial  covariance  models  (planar  and  spherical). 

The  second  part  of  his  report  dealt  with  the  following  question:  What  happens 
if  the  "true"  covariance  function  is  unknown  and  least-squares  collocation  is, 
instead,  performed  with  a "wrong"  (or  more  precisely,  non-optimal)  covariance 
function? 

How  does  the  result  change  with  respect  to  the  optimal  case,  and  what  is  the 
effect  on  accuracy  studies  ? 

Moritz  treated  the  influence  of  covariances  on  the  results  of  collocation.  For- 
mulas were  developed  for  the  standard  error  of  collocation  results  when  using  non- 
optimal  covariance  functions,  also  for  the  case  of  stepwise  collocation.  Finally 
the  behavior  of  interpolation  errors  with  and  without  the  additional  use  of  horizontal 
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gradients  was  studied  by  means  of  power  series  expansions  for  covariance 
functions  and  by  means  of  Gaussian  covariance  functions.  Moritz  concluded 
that  non- optimal  covariance  functions  had  relatively  little  influence  on  the 
interpolated  values  but  a very  strong  effect  on  covariances  as  calculated  using 
the  conventional  formulas. 

— 3 Least-Squares  Collocation  as  a Gravitational  Inverse  Problem 

In  geodesy  and  geophysics  we  frequently  meet  with  the  situation  that  a 
model  defined  by  a set  of,  say,  N parameters  is  to  be  determined  from  a 
smaller  number  n N of  observations. 

As  an  example,  the  internal  structure  of  the  earth  may  be  defined  by  a set 
of  N Parameters  describing  the  density,  the  rigidity,  and  the  compressibility 
of  the  earth  as  a function  of  depth.  The  n observations  comprise  velocities  of 
seismic  surface  waves,  together  with  the  mass  and  the  polar  moment  of  inertia 

of  the  earth.  If  the  model  for  the  earth's  internal  structure  is  to  be  realistic 
then  N will  be  large  and  n < N. 

We  thus  have  n N equations  for  N unknowns,  which  is  obviously  an 
underdetermmed  problem  admitting  an  infinite  number  of  possible  improperly 

posed  problem.  (A  problem  is  properly  posed  if  it  has  a unique  solution  that 
depends  continuously  on  the  data. ) 

Originally,  the  equations  expressing  the  data  x,  as  functions  of  the  model 
parameters  sr  will,  in  general,  be  nonlinear; 

x f - f , ( Sj , s2 , . . . , Sn  ),  i = l,2,...,n. 

By  a suitable  application  ot  Taylor’s  theorem  it  is  usually  possible  to  approximate 
these  equations  by  linear  ones: 

N 

^ 11  1 r s r » 


X , = 


or  in  matrix  notation: 


r =t\ 


x = A_s  . 

The  formal  solution  of  this  system  of  linear  equations  may  be  written  as 

- = A~ 1 * ■ 

“ ~ "Tre  a rc'g“'"  s<|uart'  matrlx'  th«  A'  1 would  be  the  ordinary  inverse  matrix 

of  A.  In  our  underdetermined  case,  however,  A'1  must  be  understood  in  the 

sense  of  generalized  matrix  inverses. 


At  any  rate,  the  solution  of  the  above  equations  may  be  considered  as  an 
inversion  of  these  equations  with  respect  to  the  parameters  sr  , which  accounts 
for  the  name,  geophysical  inverse  problems. 

Another  typical  example  of  an  "improperly  posed"  inverse  problem  is  the 
determination  of  subsurface  mass  distributions  which  produce  a given  anomalous 
gravity  field  at  the  earth’s  surface.  This  problem  is  sometimes  called  an  inverse 
problem  of  potential  theory. 

The  determination  of  the  earth's  external  gravitational  field  from  geodetic 
gravimetric  and  satellite  data  may  also  be  considered  as  an  inverse  problem 
that  is  mathematically  quite  similar  to  the  determination  of  the  internal  structure 
of  the  earth  from  seismic  and  other  data. 

This  geodetic  inverse  problem  is  likewise  underdetermined.  The  external 
gravitational  field  requires  for  a complete  description  an  infinite  number  of 
parameters,  for  instance,  the  set  of  all  coefficients  in  the  expansion  of  the 
external  gravitational  potential  in  spherical  harmonics.  This  infinite  number, 

N = 00  , of  parameters  is  to  be  determined  from  a finite  number  n of  obser- 
vations. 

Even  in  the  seismic  inverse  problem  it  is,  at  least  theoretically,  appropriate 
to  take  N = ^ if  we  wish  to  admit  reasonable  general  functions  for  density, 
rigidity,  and  compressibility  because  it  cannot  be  assumed  a priori  that  such 
functions  depend  on  a finite  number  of  parameters  only. 

Thus,  in  general,  the  space  of  parameters  will  be  infinite -dimensional 
rather  than  N-dimensional.  In  other  words,  the  proper  general  setting  for 
(linear)  geodetic  and  geophysical  inverse  problems  will  be  infinite-dimensional 
Hilbert  space.  This  was  pointed  out  by  Krarup  (1969)  for  the  geodetic  case  and 
by  Backus  (1970)  for  geophysical  inverse  problems. 

The  geodetic  inverse  problem,  the  determination  of  the  external  gravita- 
tional field  from  data  of  different  kind,  is  usually  solved  by  least-squares  col- 
location. This  technique  has  many  features  in  common  with  other  geophysical 
inversion  methods.  It  may,  therefore,  be  of  interest  to  compare  these  techniques 
and  to  exhibit  some  cross-connections.  This  was  done  by  Moritz  in  his  report 
(1976b):  "Least-Squares  Collocations  as  a Gravitational  Inverse  Problem.  " 

He  also  discussed  least-squares  collocation  from  the  point  of  view  of 
analytically  representing  the  external  gravitational  field  by  a linear  combination 
of  suitable  simpler  harmonic  functions. 

Moritz  said  that  the  subject  of  his  report  is  purely  conceptual,  aiming  at 
a better  understanding  of  least-squares  collocation  by  considering  it  in  relation 
to  other  methods;  no  new  computational  formulas  were  derived.  His  paper  is 
useful  as  a contribution  to  the  present  discussion  on  the  conceptual  foundations  of 
least-squares  collocation. 
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1.4  On  the  Computation  of  a Global  Covariance  Model 


In  his  report  Moritz  ( 1976a)  showed  that  the  local  behavior  of  a covariance 
function  of  the  gravity  anomaly,  Ag,  can  be  characterized  quite  well  by  means 
of  three  constants:  C0  variance  of  £g,  £ correlation  length  and  G0 
variance  of  the  horizontal  gradient  of  Ag.  In  his  rc|K>rt  "On  the  Computation 
of  a Global  Covariance  Model"  (Moritz,  1977),  he  treated  the  problem  of  de- 
termining a global  covariance  function  if  the  following  data  are  given:  the  var- 
iance of  gravity  anomalies,  the  variance  of  second-order  gradients,  the  corre- 
lation length  and  the  lower  degree  variances.  He  proposed  that  covariance  model 
is  a linear  combination  of  the  reciprocal  distance  covariance  function  and  a co- 
variance  function  of  logarithmic  type.  He  also  found  tint  the  correlation  length 
is  already  fixed,  within  rather  narrow  limits,  by  the  remaining  data.  It  turned 
out  that  the  gradient  variance,  G0,  rather  than  the  correlation  length,  £, 
provides  a basis  for  calculating  a global  covariance  function  from  local  covari- 
ance functions. 


1.5  Statistical  foundations  of  Collocation 

Users  of  least-squares  collocation  ask  for  a theory  that  gives  an  answer  to 
practically  meaningful  questions:  What  is  the  accuracy  of  our  results  ? Can  we 
apply  statistical  testing  techniques?  How  can  w'e  compute  statistical  distributions 
of  gravity  anomalies  or  of  deflections  of  the  vertical?  A reasonable  answer  to 
these  questions  requires  some  statistical  theory  of  the  anomalous  gravitational 
field.  But  is  this  field  really  a stochastic  phenomenon?  Such  questions  seem  to 
motivate  research  into  the  statistical  foundations  of  collocation. 

Least-squares  collocation  has  its  roots  in  many  fields; 

1.  Least-squares  estimation; 

2.  Prediction  theory  of  stochastic  processes; 

3.  Approximation  theory; 

4.  Functional  analysis,  especially  the  theory  of  Hilbert 
spaces  with  kernel  functions; 

5.  Potential  theory; 

6.  Inverse  and  improperly  posed  problem. 

All  of  these  "many  facets  of  collocation"  present  relevant  aspects  which  must  be 
taken  into  account  in  a complete  and  balances  treatment. 

The  relation  to  the  theory  of  inverse  problems  is  clear;  our  data  are  function- 
ally related  to  the  gravitational  field;  to  determine  this  field  from  the  data,  we 
must  somehow  invert  those  functional  relations.  Now’  the  gravity  field  requires 
infinitely  many  parameters  for  its  full  determination;  the  number  of  measurements, 
however,  is  essentially  finite.  Therefore,  we  have  an  improperly  posed  problem. 

To  get  a unique  solution,  we  must  impose  additional  conditions,  which  may  have 
the  form  of  a least-squares  principle  or  of  a norm  in  Hilbert  space. 
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Historically,  collocation  has  developed  from  least-squares  prediction  of 
gravity  anomalies,  which  is  an  application  of  the  prediction  theory  of  stochastic 
processes.  Hence,  statistical  considerations  l ave  played  an  essential  role  in 
collocation  from  the  ve-y  beginning. 

Also,  the  relation  to  classical  least-squares  adjustment  has  soon  been 
noted.  In  fact,  collocation  models  bear  formal  resemblance  to  conventional 
adjustment  models.  The  characteristic  difference,  however,  is  the  infinite 
number  of  parameters  necessary  to  fully  characterize  the  gravitational  field. 

This  fact  furnished  an  essential  link  to  stochastic  processes  and  to  infinite  - 
dimensional  Hilbert  spaces. 

Least-squares  estimation  and  stochastic  processes  give  a very  convenient 
mathematical  formalism  and  terminology.  They  also  provide  the  basis  fora 
statistical  interpretation  of  the  results,  essential  for  feasibility  studies. 

The  practical  success  of  the  statistical  treatment  of  collocation  has  some- 
times overshadowed  its  equally  significant  analytical  aspects,  especially  the  fact 
that  there  is  a clean  analytical  structure  underlying  it.  This  mathematical 
structure  Is  based  on  the  harmonic  character  of  the  anomalous  gravitational 
field  and  on  the  fact  that  all  quantities  of  this  field  can  be  expressed  as  linear 
functionals  of  the  anomalous  potential.  The  analytical  character  of  collocation 
is  best  brought  out  by  approaching  it  from  the  standpoint  of  approximation  theory, 
working  in  a Hilbert  space  with  a kernel  function. 

These  two  aspects,  the  statistical  and  the  analytical  aspect,  are  both  indis- 
pensable and  mutually  complement  each  other.  In  fact,  evident  already  in  the 
fundamental  paper  (Krarup,  1969),  seems  to  be  generally  agreed  upon,  although 
there  is  some  controversy  on  details. 

A literal  interpretation  of  the  anomalous  gravitational  field  as  a stochastic 
process  has  encountered  two  objections.  First,  there  is  only  one  Earth;  a 
probability  space  of  many  possible  earths  is  logically  unobjectable,  but  appears 
unnatural,  since  all  realizations  except  one  (the  real  Earth)  are  unobservable. 

Secondly,  Lauritzen  (1979)  has  proved  that  there  is  no  ergodic  Gaussian  process, 
harmonic  outside  a sphere.  This  has  sometimes  been  misinterpreted  as  a proof 
that  no  ergodic  process  modeling  the  anomalous  gravity  field  exists  at  all,  so 
that  the  covariance  function,  in  principle,  cannot  be  estimated  from  the  data. 

In  fact,  however,  the  Gaussian  structure  enters  essentially  into  Lauritzen's 
proof,  and  there  do  exist  non-Gaussian  ergodic  processes  suitable  for  collocation. 

In  his  paper  "Statistical  Foundation  of  Collocation"  (1978a),  Moritz  deals 
with  mathematical  models  suitable  as  a basis  for  statistical  treatment  of  collo- 
cation. As  a preparation  he  discusses  first  the  stochastic  processes  on  the 
circle,  because  such  processes  arc  simple  to  understand  and  they  exhibit  already 
essential  features  of  the  problem..  Then  he  treats  stochastic  processes  on  the 
sphere,  which  are  suitable  as  statistical  models  for  collocation.  Moritz  shows 
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that  Lauritzen's  theorem  on  the  non  ex  is  lance  of  orgodic  Gaussian  stochastic 
process  models  for  collocation  to  be  essentially  dependent  on  the  Gaussian 
character.  lie  gives  two  non-Gauss ian  crgodic  models,  one  of  a ginuinel\ 
probabilistic  character  similar  to  1 aurit/.en's  model,  and  another  based  on  a 
formal  probability  theory  in  rotation  group  space. 

1 his  second  model  gives  a ta'istieal  foundation  of  the  usual  homogeneous 
and  isotropic  covariance  anal\s  1 tht  nomalous  gravity  field.  The  model  also 
provides  a basis  for  the  study  of  th<  ■ ii  'lie  d distribution  of  quantities  related 
to  this  field.  It  also  allows  a form:  ! s*.m:stical  treatment  of  the  anomalous 
gravitational  field  which  is  indepei'  : 1 m interpretation  of  this  field  as  some 

genuinely  physical  stochastic  pro.-*  . < iv 5 seems,  therefore,  to  be  preferable. 

If  the  approach  presented  by  ' >r  is  accepted,  then  a detailed  theory  of 
statistical  distributions  for  ge<>d<  I!y  want  quantities,  such  as  gravity 
anomalies,  geoidal  heights,  and  feet:  ms  of  the  vertical,  could  be  developed 
and  applied  to  the  statistical  testing  " results  of  least-squares  collocation. 


1.6  The  Operational  Approach  to  ! h . sical  Geodesy 

There  are  essentially  two  pos  tbit  approaches  to  physical  geodesy  (as  also 
to  other  natural  sciences):  they  might  be'  called  the  model  approach  and  U.e 
operational  approach.  Essentially,  the  first  approach  starts  from  a theory  , the 
second  from  the  observations.  Obviously,  the  two  approaches  are  closely  related 
to  the  deductive  method  and  the  inductive  method  in  the  natural  sciences. 

In  the  model  approach,  one  starts  from  a mathematical  model  or  from  a 
theory  and  then  tries  ,o  fit  this  model  to  re  dity,  for  instance  by  determining  the 
parameters  of  this  model  from  obse rvapnns.  The  classical  geodetic  example  are 
the  centuries-old  attempts  to  determine  the  parameters  of  an  earth  ellipsoid  In 
observations. 

Perhaps  the  most  elaborate  form  of  this  model  approach  is  the  boundary- 
value  problem  of  physical  geodesy  in  (he  formulation  of  Molodenskv.  It  has  a 
mathematically  enormously  interesting  an  i deep  theory  and  is  practically  highly 
significant,  as  the  many  gravimetric  gcoid  determinations  and  computations  of 
deflections  of  the  vertical  shov*.  Howvvi  t . this  approach  has  ils  weaknesses:  the 
required  continuous  gravity  coverage  is  practically  not  realizable;  on  the  other 
hand,  many  other  important  data  cannot  be  incorporated  into  this  theory.  The 
model  selects  its  data. 

At  present  we  have  a great  number  of  geodetic  measurements  of  very 
different  types,  from  terrestrial  angle  and  distance  measurements  to  satellite 
data  of  various  kinds.  The  question  arises;  how  can  we  use  and  combine  all 
these  data  in  the  best  possible  way  Phis  iv  the  operational  approach. 


J 


Let  us  summarize.  In  the  model  approach  one  asks:  how  can  I best  de- 
termine my  model  by  suitable  observations?  In  the  operational  approach  one 
asks;  how  can  I make  best  use  of  all  my  observations? 

As  a matter  of  fact,  the  two  approaches  do  not  compete  with  each  other; 
each  one  incorporates  important  aspects,  and  the  two  approaches  mutually 
complement  each  other. 

l’he  operational  approach  to  phy  sical  geodesy  has  come  up  at  a relatively 
recent  date,  when  a huge  number  of  measurements  of  new  types  was  available 
and  when  it  turned  out  that  the  classical,  especially  the  gravimetric  approach 
failed  to  give  a complete  answer  in  view  of  the  lack  in  gravity  data. 

In  geometrical  geodesy  already  least-squares  adjustment  is  in  the  spirit 
of  an  operational  approach  (how  can  1 best  use  all  my  measurements).  In 
physical  geodesy,  operational  methods  have  been  known  under  the  names  "least- 
squares  collocation,  " "integrated  geodesy,  " "operational  geodesy'."  All  these 
methods  are  very  similar;  they  all  aim  at  an  adequate  treatment  of  the  gravity 
field,  in  addition  lo  an  adjustment  of  measuring  errors.  They  all  use  quadratic 
minimum  principles  incorporating  not  only  the  measuring  errors,  but  also  the 
anomalous  gravity  field. 

Moritz  (1978b)  presented  a systematic  treatment  of  the  operational  approach 
in  his  report:  "The  Operational  Approach  to  Physical  Geodesy."  He  shows  that, 
using  contemporary  mathematical  techniques,  a straight  road  leads  from  the  non- 
linear observational  equations  to  collocation  with  kernel  functions  and  least-squares 
collocation.  After  linearization  an  improperly  posed  problem  was  obtained,  to 
which  Moritz  applied  three  different  standard  methods  of  solution:  1 ) a restriction 
of  the  solution  space,  leading  to  "pu  re -collocation,  " 2 ) variational  principles 
of  Tichonov  type,  by  which  measuring  errors  can  be  taken  into  account,  leading  to 
a generalized  collocation  w ith  kernel  functions,  and  3 ) a statistical  approach, 
leading  to  least-squares  collocation. 

At  the  end  of  his  report,  Moritz  discussed  these  various  stages  and  possible 
alternatives.  All  ihroe  approaches  seem  to  converge  on  collocation  with  kernel 
functions  and  least-squares  collocation. 

1 •1 * * * * *  7 A ppr< uomation  of  Covariance  Functions  by  Non-positive  Definite  Functions 

During  the  last  decade,  when  least-squares  collocation  presented  itself  as 

the  data  processing  model  in  physical  geodesy,  the  most  serious  argument  against 

was  the  inversion  of  a large  matrix  resulting  in  much  computer  time  needed  for 

this  purpose.  No  attention  was  paid  to  the  time  used  for  calculating  the  linear 

functionals  on  the  covariance  function  because  for  simple  pi'oblems  this  time  is 
definitely  inferior  to  the  inversion  time.  The  situation,  however,  changed  im- 

mediately when  problems  were  attacked  which  involved  many  and/or  difficult 
covariance  calculations.  Although  Itapp  and  Tsc homing  (1974)  succeeded  in 
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deriving  closed  expressions  of  covariance  functions  for  different  models  of 
anomaly  degree  variances,  the  closed  expressions  still  consist  of  functions  like 
logarithmic  and  trigonometric  functions,  which  are  expensive  in  terms  of  com- 
puter time. 

We  mention  only  a few  kinds  of  application:  prediction  of  mean  gravity 
anomalies  over  rectangular  blocks  from  point  values,  prediction  of  mean  gravity 
anomalies  over  larger  areas  from  mean  gravity  anomalies  over  smaller  areas, 
prediction  of  mean  gravity  anomalies  from  satellite  altimetry  data,  all  problems 
involving  satellite  dynamics.  All  these  applications  have  one  common  feature: 
it  is  necessary  to  calculate  covariances  by  numerical  integration.  In  case  of 
mean  gravity  anomaly  prediction  the  integration  is  at  most  twofold,  in  case 
of  satellite  dynamics,  however,  it  is  multifold,  in  the  former  case  an  explicit 
integration  procedure  can  be  avoided,  if  one  replaces  the  rectangular  area  of 
integration  by  a circular  one.  The  so-called  smoothing  opt* ration  is  caused  by 
an  isotropic  smoothing  operator  acting  on  the  covariance  function  which  itself 
is  also  isotropic.  Therefore,  the  convolution  of  the  smoothing  operator  with 
the  covariance  function  corresponds  to  a product  of  the  corresponding  eigen- 
values, which  is  naturally  very  simple.  In  order  to  obtain  closed  expressions 
for  the  mean  gravity  anomaly  function,  however,  a further  artificial  assumption 
has  to  be  made:  the  eigenvalues  of  the  smoothing  operator  have  to  be  replaced 
by  some  other  values  (Schwarz,  1976).  In  problems  involving  satellite  dynamics 
as  satellite-to-satellite  Kinging  probably  the  only  way  to  calculate  covariances 
is  by  numerical  integration  over  some  time  interval.  Using  exact  covariances 
for  the  integration  procedure  is  extremely  time  consuming. 

These  were  the  reasons  why  the  question  arose  whether  it  is  possible  to 
use  some  more  or  less  accurate'  approximations  of  the  exact  covariance  functions; 
the  approximating  functions  should  be  simple,  easy  to  handle,  accurate  and  should 
consume  as  little  mass  storage  as  possible.  Sunk  el  (1978)  studied  this  problem 
and  reported  his  findings  in  the  report;  "Approximation  of  Covariance  Functions 
by  Non-Positive  Definite  Functions.  " He  studied  three  different  kinds  of  approx- 
imating functions,  all  of  them  being  finite  elements;  1 ) a step  function,  2 ) 
a piecewise  linear  function  and  3 ) a cubic-spline  function. 

The  basic  principle  underlying  his  investigations  was  well  known  and  fre- 
quently applied  in  many  fields;  the  network  principle;  generate  a net  of  fixed 
points  (here  grid  points)  and  perform  very  accurate  measurements  at  these 
points  (here,  calculate  exact  covariances);  these  fixed  points  serve  as  a basis 
for  small  scale  measurements  which  can  be  performed  using  simpler  apparatus 
(here,  more  or  less  interpolation  of  covariances  by  means  of  finite  elements). 

His  report  was  primarily  devoted  to  the  study  of  interpolation  errors,  pertur- 
bation of  spectra  and  to  the  consequences  of  the  approximation  for  the  predicted 
signal  and  its  mean  square  error.  Sunkel  concluded  that  because  of  its  smooth- 
ness and  its  most  favourable  approximation  properties  the  spline  function  rep- 
resentation of  the  covariance  function  presents  itself  as  a very  useful  tool  for  this 
kind  of  application. 
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1.8  Covariance  Expressions  for  Second  and  Lower  Order  Derivatives  of  the 


Anomalous  Potential  ' 

Since  Moritz  (1972)  discussed  the  theory  and  application  of  least  squares 
collocation  techniques  for  gravimetric  geodesy  problems  a number  of  additional 
studies  have  been  done  to  gain  better  insight  into  collocation  from  a theoretical 
point  of  view,  and  to  apply  least  squares  collocation  techniques  to  geodetic  prob- 
lems of  interest. 

One  aspect  of  this  work  is  the  development  of  consistent  covariance  func- 
Uons.  One  step  in  this  direction  was  described  in  Tscherning  and  Rapp  (1974). 

In  this  paper  an  anomaly  degree  variance  model  was  developed  that  was  fitted 
to  different  types  of  information  that  could  be  related  to  the  anomaly  degree 
variances.  In  that  paper  the  following  model  was  chosen: 


l-  ( a - 2 ) ( i + b ) ' ' 

where  A = 425.28  mgal2,  B - 24,  and  a parameter  s (used  in  the  application 
of  Cj,  ),  s = 0.999617.  In  that  report  a number  of  covariance  and  cross  cov- 
ariance functions  were  derived  involving  anomalies,  deflections  of  the  vertical, 
and  geoid  undulations.  All  these  quanties  were  considered  to  be  global  covariance 
functions  based  on  a single  anomaly  degree  variance  model. 

What  was  lacking  from  the  above  study  were  the  covariances  involving  a 
number  of  derivatives  that  could  be  related  to  gravity  gradients.  To  extend  the 
earlier  work  Tscherning  (1976)  developed  an  extended  set  of  covariance  function 


equations  for  use  with  several  different  anomaly  degree  variance  models.  More 
specifically  Tscherning  (ibid.)  chose  to  model  the  potential  degree  variances 
which  refer  to  the  disturbing  potential.  Thus  if  cov(TP,TQ)  is  the  covariance 
between  the  disturbing  potential  at  points  P and  Q he  writes: 


cov  ( Tp , Tq 


> -y  • 


c 0 + i 

£<T,T)S*  P. 


( T,T) s 


where;  a ^ (T,T)  are  the  model  potential  degree  variances  and  cr^  (T,T) 
are  corrections  to  the  model  potential  degree  variances  to  degree  n.  We  have; 


s R„7(rp*rQ) 


t = COS  *1 > 


where  Rb  is  the  radius  of  the  Bjerhammer  sphere,  r is  a geocentric  distance 
and  i b is  the  arc  between  the  two  points.  Tscherning  (ibid.)  writes  the  O (T,T) 
values  in  the  following  general  form: 

Of  = A,  n ( t + kj)-1  (4) 


J =0 
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where  i (1,  2,  or  3)  is  a model  number  and  A.  is  a constant  in  units  of 
(m/see)'1 . The  a;  values  are  related  to  the  e values  In  (Tscherning  and 
Rapp,  ibid. ) 

h.: 

a?  = jr^rr  c>  (5> 


Using  (4)  and  (5)  with  (1)  Tschcrnm  gave;  k 
A = A,  10lo/H„8  with  It t,  in  meters. 


-2,  kj  1 , kc 


At  this  point  Tschcming  (197b)  developed  tin1  necessary  equations  to  eval- 
uate the  covariances  and  auto  eovariani  '•  r the  following  qunntitics;(  ] ) the 
height  anomaly;  (2)  the  negative  graviu  < suirbanee  r • (!i)  the  gravit\ 
anomaly  ( Ag) ; ( 4)  the  radial  compo  •>  nt  of  the  gradient  of  A ; ; ( r»)  the 
second  order  radial  derivative  >f  ; (i  i,(7)  the  latitude  and  longitude  com- 
ponents of  the  deflection  of  th< • vert  ■!-,  i s),r<»)  the  derivatives  in  the  northern 
and  eastern  direction  of  Ag;  ( 10  » , ( 1 1 the  derivatives  of  the  gravity  distur- 
bance in  the  northern  and  eastern  din  ction;  ( 12)-(  14  ) the  second  order  de- 
rivatives of  T in  the  northern,  in  mixed  northern  and  eastern,  and  in  the  eastern 
direction.  Specific  equations  for  the  , giant ities  can  lx  found  in  I'schcming 
(ibid. , p.  2 and  3). 

To  implement  these  equation.'-  rsche ruing  devised  a subroutine  called 
COVAX  that  was  given  in  the  report.  A few  corrections  to  the  text  and  the 
computer  program  have  been  made  since  its  original  publication.  Such  correc- 
tions can  be  obtained  directly  from  Tscherning. 

1.9  Two  Models  for  the  Degree  Variain-  of  Global  Covariance  Functions 

In  carrying  out  the  computations  wth  COVAX,  Tscherning  found  that  the 
variance  of  the  vertical  gradient  of  the  gravity  anomaly  was  about  70001  " 
which  implies  a horizontal  anomaly  gradient  (C  ) variance  of  3500  K . Moritz 
(1977)  showed  that  this  high  variance  implies  a correlation  length  (()  smaller 
than  found  in  practice.  He  then  postulated  a form  of  an  anomaly  covariance 
function  that  implied  an  anomaly  degree  variance  model  that  would  avoid  the 
apparent  problems  of  the  Tscherning  Rapp  model  when  dealing  with  gradient 
covariance1  functions.  Jekeli  ( 3 9 7 s ) proceeded  to  investigate  the  Moritz  sug- 
gestion and  attempted  to  provide  new  numerical  estimates  of  iht  anomaly  degree 
variance  models. 


Jekeli  (ibid.)  expressed  the  anomaly  covariance  function,  (’(  in 

the  following  form:  ^ 

C(I\ti)  = ^ e,  (— ^ — P (cos  U)  (6) 

£ OP  r.;  / l 


where  Rc  is  the  radius  of  a sphere  to  which  the  e values  are  referred  and 
is  given  by:  . . 


cz  = n 1 (T+T)cr^  + *“77^2 


( f - 2)(  P + 14) 


a^4  ’,  f 3 


In  this  expression  the  parameters  to  be  estimated  are  or,  , oa  , crl  , and  cr2  . 

It  was  noted  that  i)  ce,  0,  e(|nation  (7)  reduces  to  the  Tscherning-Rapp  model. 
The  values  of  A and  li  are  not  easils  found  through  adjustment  procedures; 
rather  a trial  and  error  technique  is  chosen. 

■lekeli  then  proposed  to  determine  the  parameters  of  the  model  by  using 
the  following  data;  anomaly  degree  variances  to  degree  52  from  Rapp  (1977)  or 
to  degree  20  from  the  (iKM  > solution;  point  anomah  variance;  the  vertical 
gradient  variance;  and  mean  anomaly  variances  for  1°  and  5°  blocks.  Special 
consideration  was  given  to  using  a high  gradient  variance  (7000 Es)  and  a low 
gradient  variance  (4001"')  to  see  the  effect  on  the  adjusted  model. 

Many  different  solutions  were  made  by  .lekeli  to  obtain  solutions  for  the 
two  eomponi  at  (o,  , cy  a 0)  model  and  a one  component  (a-,  = 0,  0/ a J 0 ) model. 
Using  the  i,i-  \]  9 degree  variances  to  degree  20  the  results  of  Jekeli’s  model 
fits  arc  shown  in  Table  1 (from  .lekeli.  p.  53).  The  best  fit  to  the  data  is  found 
to  be  with  the  two  component  model  with  the  low  gradient  variance.  The  one 
com|x>nent  model  does  not  give  a good  (joint  anomaly  variance  when  a low  gra- 
dient varance  is  used.  At  the  lower  degree  (?  20)  the  c,  value  from  one 

model  can  be  twice  that  of  another  model.  At  degree  60  this  difference  has 
decreased  significantly , the  magnitude  of  c,  being  on  the  order  of  5 mgal^  for 
all  models.  These  values  are  somewhat  high  when  compared  to  recent  (Rapp 
197H)  determinations  showing  values  more  on  the  order  of  3 mgal2. 

1-10  A Model  Comparison  in  Least-Squares  Collocation 

In  addition  to  considering  improved  anomaly  degree  variance  models  needed 
for  covariance  function  computation,  several  studies  were  carried  out  to  obtain 
a better  understanding  of  some  theoretical  problems  in  least  squares  collocation, 
one  such  study  was  carried  out  by  Hummel  (1976).  Here  Rummel  pointed  out 
tliat  two  models  of  least  squares  collocation  arc  in  current  use.  The  first  he 
refers  to  as  the  model  of  Method  One  which  is: 

l A x + s ’ + n (8) 

where;  l is  a vector  of  observations; 

x is  a vector  of  unknown  parameters; 

A is  the  coefficient  matrix; 

s’  is  the  random  signal  part  of  l ; 

n is  a random  signal  or  noise. 

Ih<  estimation  of  the  signal  at  an  arbitrary  point  is  then  given  by  (Rummel,  ibid) 

■s  <'»«'(  r,','  + Cnn)  (2-  A x)  (9) 

vvh‘,|V  r **  is  the  covariance  matrix  between  the  "observed”  quantities  and 
c“'  is  the  e (Variance  matrix  between  the  signal  (s)  to  be  estimated  and  the 
’observed"  signal.  Model  One  is  the  model  originally  suggested  by  Moritz  (1972). 


Table  1.  Parameters  of  Anomaly  Degree  Variance  Models  (from  Jekeli,  1978). 


>1  211:  two-component  model  (3.23)  with  an  observed  value  of  Go  h **  3500  E3 
(1  1L:  single  component  model  (3.12)  with  an  observed  value  of  G0h  ■=  200  E 
1 lHt  single  component  model  (3. 12)  with  an  observed  value  of  G0  m “ 3500  E' 


A second  model  for  use  in  least  squares  collocation  was  suggested  by 
Merit  ’ and  Schwarz  (197,'}).  This  model  is  called  Model  Two  by  Hummel  and  is 
written  in  the  following  form: 


t Ax+lts  + n (10) 

where  s a the  desired  signal  and  appears  directly  in  the  model.  The  values 
of  s,  x,  md  n in  (10)  are  not  necessarily  the  same  as  the  values  of  in  (8)  and  (9). 
The  solution  for  a signal  is  given  by  (Hummel); 

s C„  R'(HC„RT +Cn.)"1  (*- Ax)  (11) 

It  is  this  form  of  the  collocation  solution  that  is  often  reduced  to  one  involving 
an  inversion  of  a matrix  whose  size  is  equal  to  that  of  the  signals  being  estimated, 
and  not  on<  requi ring  the  inversion  of  a matrix  whose  size  corresponds  to  the 
number  of  observations  used  in  the  process. 

Hummel  clearly  shows  that,  in  general,  the  results  from  (9)  and  (11)  will 
not  be  the  same.  The  onh  case  of  equivalence  will  be  when  we  have; 


•s'  = H s (12) 

In  most  gravimetric  applications  H will  have  to  be  a matrix  transforming  an 
infinite  set  of  signals  s (such  as  gravity  anomalies)  into  s'  (such  as  a geoid 
undulation).  In  practice  this  could  only  be  done  in  an  approximate  way.  Thus 
Rumrnel  proved  that  the  two  models  used  in  least  squares  collocation  are  differ- 
ent. Only  in  special  cases  should  Model  Two  be  preferred  with  most  applications 
being  done  using  Model  One. 

1.11  Potential  Coefficient  Determinations  from  10c  Terrestrial  Gravity  Data 
By  Means  of  Collocation 


A specific  application  of  Model  One  was  carried  out  by  Sjoberg  (1978a) . 
tie  considered  the  estimation  of  potential  coefficients  from  mean  gravity  anomaly 
data  given  in  41fi  10°  equal  area  anomaly  blocks.  In  this  case  equation  (9)  takes 
the  following  fornit 


if;}  * 


(13) 


where  C,  S are  the  predicted  fully  normalized  potential  coefficients;  C is  they 
covariance  matrix  between  the  mean  anomalies  and  cCfcs  are  the  cross  covari-  y 
ances  between  the  potential  coefficients  and  the  mean  anomalies;  D is  the  e^ror.  jr 
covariance  matrix  of  the  anomalies.  Instead  of  computing  the  mean  covariano^4^  / 
functions  by  numerical  integration  Sjoberg  (ibid)  developed  several  diffident  f 

procedures  including  the  use  of  the  smoothing  operator  of  Pellinen,  anJrthe 
evaluation  of  the  point  covariance  functions  at  a certain  elevation  above  the  mean 
sphere. 
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Sjoberg  carried  out  three  diffe rent  solutions:  1 ) a least  squa  rc“s  eollocat ion 
solution  considering  the  mean  elevations  of  (he  anomaly  blocks;  2 ) a least 
squares  collocation  solution  setting  the  elevations  to  zero;  and  3 ) a collocation 
solution  setting  the  elevations  to  x<  m and  the  noise  matrix  1)  to  zero.  These 
solutions  were  then  compared  to  the  potential  coefficients  found  from  the  usual 
numerical  integration  proeedu re,  t. . the  OKM  <)  potential  coefficients,  and  to  co- 
efficients computed  from  5'  mean  tv  malies  (Rapp,  1977a).  Some  re  tills  in 
terms  of  differences  are  summ  tn  ' in  table  2. 


Table  2. 

Comparison  of  If 

•!  n!  I ( Oeffieients 

From  Least 

Squares  Co 

1 location  and  < ■;  In 

al  ■ lata  from  S j 

joberg,  1978a) 

Pereen'ag- 

MS  Fndulation  1 

i RMS  Anomalv 

| Difforoiu  ( fb  once  (m) 

Difference  (mgals) 

1 

GEM  9,  1*  251 

1 

s.  7 

7.2 

GEM  9,  Int.  10° 

73 

8.  8 

7.6 

GEM  9,  Coll.  1 

68 

<),  2 

7.  2 

GEM  9,  Coll.  2 

GH 

x . 8 

7.2 

GEM  9,  Coll.  3 

71 

9.9 

7.5 

P 251,  Int.  10 

32 

2.  0 

3. 6 

# 251,  Coll.  1 

30 

2.  8 

3.  3 

# 251,  Coll.  2 j 

30 

3.0 

3.3 

Int.  10°,  Coll.  1 ' 

13 

2.  1 

1.6 

Int.  10°,  Coll.  2 1 

14 

2.  6 

1.7 

Int.  10°,  Coll.  3| 

15 

3 2 

1.7 

Coll.  1,  Coll.  2 ! 

3 

2. 1 

0.7 

From  this  data  several  things  can  bo  s,  , n. 

1.  The  coefficients  from  the  r,  solution  agree  better  with  the  (JEM  9 coeffi- 
cients than  the  10°  solutions. 

2.  The  effect  of  including  fL  ••  ms  in  the  tin  collocation  solution  slightly 

degrades  the  comparisons  with  Of  M 

3.  The  best  agreement  with  the  coefficients  from  the  10°  anomalies  hikes 
place  when  the  D matrix  is  included  with  the  collocation  solution. 

4.  The  collocation  solutions  agree  better  with  the  GEM  9 coefficients  than 
do  the  coefficients  found  from  the  usual  integration  procedures. 

Sjoberg  (ibid)  also  considered  the  computational  effort  in  carrying  out  these 
computations.  He  found  that  19  seconds  of  computer  time  were  needed  for  the 
usual  integration  (or  summation)  procedure  while  561  sebonds  were  needed  for  the 
collocation  2 solution.  Although  the  coefficients  from  the  collocation  solution  may 
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1)0  slightly  hotter  than  tho  summation  the  tremendous  increase  in  computational 
effort  does  not  justify  the  use  of  least  squares  collocation  for  the  estimation  of 
potential  coefficients  from  a global  set  of  gravity  anomalies. 

1.  1-  A Comparison  of  Bjcrha  mmar's  Method  and  Collocation  in  Physical  Geodesy 

Another  study  was  conducted  by  Sjoberg  (1978b)  on  the  relationship  between 
least  squares  collocation  and  the  Bjerhammar  theory.  The  Bjerhammar  theory 
postulates  a set  of  fictitious  gra\  it\  anomalies,  Ag*  , located  on  the  surface  of 
the  Bjerhammar  sphere  (which  is  internal  to  all  the  masses).  These  anomalies 
are  related  to  the  observed  anomalies  (Ag)  through  the  Poisson  integral.  In 
principal  the  Ag*  values  ean  be  found  from  solving: 

Ag  A Ag*  (14) 

where  the  elements  of  A are  determined  from  the  elements  of  the  Poisson  integral. 
The  solution  of  (1 1)  depends  on  the  number  of  given  anomalies  and  the  number  of 
anomalies  to  be  found.  A similar  equation  to  (14)  for  anomaly  prediction  can  be 
written  using  least  squares  collocation  procedures; 

Ag.  C,  (C  + l)f‘  Ag  (15) 

when  the  elements  o|  (’,  and  ( depend  on  a spatial  covariance  function.  Sjoberg 
then  examines  various  solutions  of  (14)  writing  a general  form  as: 

Ag*  = QAr(AQArf  Ag  (16) 

where  tj  is  a weight  matrix.  If  the  number  N of  the  Ag*  values  approaches 
infinity  Sjoberg  shows  that  in  this  case  the  solution  of  (16)  and  the'  resultant  pre- 
diction of  new  anomalies  from  the  original  set  can  be  written  as: 

Ag:  (’.  C _I  Ag  (17) 

where  the  elements  of  the  ('  matrix  arc  computed  from: 

r ' V ,‘‘"*'7777)  P»  (COS</'I-’)  (ls) 


Il(  re  tlie  c • are  values  uniquely  associated  w ith  a minimum  norm  process  of  a 
■4*  nf  ; ^ 1 '•  bj*  I'ba 'iim.i r method.  I’lu  s<  c ‘ values  are  analogous  to  the  anomaly 
degree  variances  (<  ) of  the  least  squares  collocation  technique.  Thus  for  each 
type  of  degree  yarianee  in  a collocation  solution  there  is  a corresponding  minimum 
norm  ' lution  In  iln  generalized  Bjerhammar  theory  (Sjoberg,  ibid,  p.  9). 

Sjoberg  (ibid)  continues  this  discussion  for  a number  of  different  cases,  such 
as  including  noise  in  the  Bjerhammar  theory,  and  an  analysis  of  a reflexive  pre- 
diotion  process  suggested  by  Bjerhammar  in  1974.  These  studies  showed  that  for 
a special  ease  of  c * en  2n  + 1 the  concept  of  reflexive  prediction  was  less 
sensitive  to  changes  4 the  radius  of  the  Bjerhammar  sphere  than  was  collocation. 
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In  the  second  method  structural  differences  between  the  gravitational  and 
inertial  fields  are  used  to  separate  the  two  effects.  The  differences  show  up  in 
the  second  and  higher  order  gradients  of  the  force  fields.  Therefore,  additional 
measurements  are  necessary  in  this  case.  Moritz  (1907)  has  shown  that  for  an 
aircraft  with  inertial  stabilization  the  second  derivatives  of  the  force  field  do  not 
contain  inertial  disturbances,  so  that  purely  gravitational  second-order  gradients 
can  be  measured.  They  are  used  to  obtain  the  gravitational  force  vector  by  in- 
tegrating along  the  flight  path.  It  should  lie  noted  that,  in  contrast  to  the  first 
method,  a rigorous  separation  of  gravitation  find  inertia  is  possible  in  this  case, 
and  that  from  a theoretical  point  of  view  this  approach  is  preferable.  The  prac- 
tical difficulties  originate  in  the  design  of  instruments  accurate  enough  to  make 
an  application  feasible.  Advances  in  instrument  development  have  been  rapid 
during  the  last  years  and  a gradiometer  with  an  accuracy  of  a few  Kotvos  may 
be  available  in  the  near  future.  Therefore,  the  capabilities  of  an  airborne 
gradiometer  system  was  studied  by  Schwarz  (197b)  in  his  report;  "Geodetic 
Accuracies  Obtainable  from  Measurements  of  First  and  Second  Order  Gravita- 
tional Graditents.  " 

The  accuracy  study  of  Schwarz  was  performed  using  the  method  of  least- 
squares  collocation.  There  were  three  reasons  why  this  approach  was  especially 
suited  for  the  problem.  First,  it  allowed  the  combination  of  heterogeneous 
data  in  a consistent  way.  This  was  very  important  because  geoidal  heights, 
gravity  anomalies,  and  different  second-order  gradients  were  used  as  measure- 
ments. They  must  be  evaluated  in  such  a way  that  their  common  origin  from  the 
same  anomalous  gravity  field  was  part  of  the  system.  In  least-squares  collocation 
this  is  achieved  by  describing  the  statistical  structure  of  the  field  by  a covariance 
function.  Second,  mean  gravity  values  at  ground  level  must  be  estimated  using 
point  values  on  a profile  in  flying  altitude  and  additional  information  on  ground. 
This  involved  interpolation  between  profiles,  downward  continuation,  combination 
of  different  quantities,  and  estimation  of  mean  values.  All  these  steps  could  be 
united  in  a single  step  procedure  in  the  collocation  method.  This  was  impos- 
sible when  using  the  corresponding  integral  formulas.  Third,  different  assump- 
tions on  the  structure  of  the  gravity  field  and  on  the  accuracy  of  the  measurements 
must  be  investigated.  Again  this  was  simple  with  the  collocation  method  because 
it  only  involved  a change  of  the  fundamental  covariance  function  or  of  the  error 
variances. 

A detailed  analysis  of  interpolation,  downward  continuation  and  mean  value 
determinations  was  given  by  Schwarz.  The  influence  of  measureing  errors  was 
considered  and  the  effects  were  discussed  in  connection  with  the  stability  prob- 
lem. Finally,  the  contribution  of  accurate  satellite  altimeter  to  a combined 
acceleromcter-gradiometer  system  was  taken  into  account. 

The  results  of  his  study  showed  that  a system  of  this  kind  could  significantly 
contribute  to  our  knowledge  of  the  anomalous  gravity  field  if  second-order  grav- 
itational gradients  could  bo  measured  with  an  accuracy  of  a few  Kotvos. 
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Sjoberg  (ibid)  found  that  the  prediction  results  and  the  stability  of  the  two  methods 
would  be  the  same  if  the  collocation  procedure  was  applied  with  half  the  depth  to 
the  Bjerhaminar  sphere  used  in  the  reflexive  prediction.  This  condition  was  based 
only  on  the  analysis  of  one  degree  varinaee  model  with  additional  analysis  being 
needed. 

2. _ Analyses  of  Airborne  Measurements 

Dr.  Schwarz  was  doing  studies  on  airborne  and  satellite  measurements  and 
their  application  for  determination  of  mean  anomalies.  He  prepared  two  reports 
and  did  extensive  computations  for  a test  area.  Lenny  Krieg  was  assisting  him  in 
the  computations  for  the  lest  area  in  the  following  Schwarz's  work  is  summarized 
based  on  his  formal  and  informal  reports.  Some  parts  of  his  reports  are  quoted 
here  without  quotation  marks. 

2.  i • . delic  Aceu nicies  Obtainable  from  Measurements  of  First  and  Second 
( bder  ( I m v itat  iona  I ( ; radients 

I >p  a imenls  w ith  airborne  gravimeters  have  been  performed  over  the  last 
17  \ car;.  M though  the  equipment  designed  for  this  objective  has  a high  degree 
of  sophisti  a lion,  the  results  obtained  so  far  are  not  accurate  enough  for  geodetic 
purposi  The  reason  lies  in  the  complicated  structure  of  the  force  field  acting 
on  i m-a.  ing  gravimeter,  (dravimeters  are  basically  accelerometers  and  they 
measu  re  the  resultant  of  gravitational  and  inertial  forces.  If  they  are  used  as 
stationary  instruments,  as  in  most  terrestrial  applications,  the  only  inertial 
force  a. -ring  on  (he  gravimeter  is  the  centrifugal  force.  Therefore,  the  output 
of  the  instrument  is  the  combined  effect  of  gravitational  attraction  and  centrifugal 
force,  be.  gravity.  The  situation  is  more  complicated  in  a moving  gravimeter. 

The  Coriolis  force  has  to  be  taken  into  account  and,  more  important,  irregular 
acceleration.;  of  the  base  w ill  strongly  influence  the  result  of  the  measurements. 
Such  undt  sir-able  inertial  forces  are  especially  strong  in  a moving  aircraft  and 
•h'  a is  i,.  .iv  to  rigorously  separate  the  gravitational  part  from  the  inertial 
pa,  i-v  sing  , ci  imeter  measurements  only.  Therefore,  additional  information 
is  necessary  to  extract  the  gravitational  effect. 

IV  in  lh<»ds  have  been  proposed  to  reach  this  goal.  In  the  first  one  infor- 
1 a ■ frequency  behaviour  of  the  different  forces  is  used  to  separate 

a* a. : ' i <i.l  disturbing  accelerations  by  statistical  filtering  techniques.  Meissl 
(IDiO.  has  investigated  this  approach  using  the  theory  of  stochastic  processes 
and  <•<  i tain  assumptions  on  the  power  spectra  of  the  force  fields  involved.  He 
concludes  that  it  is  most  difficult  to  separate  gravitation  and  inertia  in  the  me- 
dium frequency  range  with  half  wavelength  between  30  and  150  km.  This  is  only 
P'lysihb  if  detailed  information  on  the  two  spectra  is  available  which  usually  will 
not  l>e  the  ( ase.  I’he  high  frequencies  can  be  blocked  by  a low  pass  filter,  the 
lov\  frequencies  can  be  improved  by  regularly  updating  altitude  and  position.  The 
remaining  errors  will,  however,  be  of  a size  which  will  not  allow  a useful  geodetic 
application  or  the  filtered  data.  The  findings  from  a probabilistic  error  analysis 
were  confirmed  by  results  obtained  by  Szabo  and  Anthony  (1971)  in  an  analysis  of 
actual  measurements. 
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2.2  Simulation  Study  of  Ai  rheme  grid  inHri 

In  his  previous  report.  Schwarz  i I'.iT'h  studied  the  accuracy  of  airborne 
gradiometry  and  he  drew  some  <•■  n<  lus  ions  about  optimal  point  ennfigti  rot  inns 
and  data  combinations.  Ilis  rr)»ii"  Is  A • ''Simulation  Studv  of  Airborne 
Gradiometry"  supplements  some  I tv  pr<  vious  investigations.  These 
simulation  studies  displayed  t '><  r-of  individual  experiments  onh 

and  were  therefore  not  meant  ' results  of  an  accuracy  study.  I lu 

interest  of  a simulation  study  wa  ' i ■ on  of  ope  rat  ional  realization  and 
opti  ma  1 performa  nei  ■ . 


In  order  to  get  an  operation" ! < for  a irhom< ■ gradiomet  r\  the  most 

important  problem  to  eope  with  \v:  tb*  ft  ■ m m t handling  of  large  amounts  of 
data.  The  proposed  measu ring  '■ ' iroduee  about  :'.«•)  observations  p-  r 

profile  and  degree.  In  order  to  ■ ■ a ' x in  area  with  profiles  spaced  at 

1°  we  have  to  treat  120  000  me  < > nr  lor  mean  gravity  values  below 

1°  x lc  we  have  to  use  20'  spue  i i the  above  number  of  measurements 

will  triple.  It  was  shown  by  Schwa  ; : ,'7)  how  the  number  of  observations 
can  be  reduced  without  signifiean;  ;>  tiring  the  aeeu racy  of  the  results.  Fee 
an  operational  program,  however.  1 • . • --sa  ry  to  use  all  information  a \ ail 

able.  Not  so  much  to  increase  a<  eu  ta<  but  to  make  results  more  reliable. 
Then* fore,  Schwarz  incorporated  m e-w points  in  the  program. 

The  simulation  studies  presented  l Selnvarz  showed  that  least -squares 
collocation  offers  an  adequate  mode!  to  c Miniate  gravity  anomalies  from  air- 
borne grad iome ter  measurements.  Tin  pmeeduiv  is  simple  numcrieallv 
and  allows  to  handle  large  amoun"  of  ■ >th  regular  requirements  on  core 

strage  and  small  demands  on  computer  He  . 

The  deviations  of  the  estin  o d vitv  anomalies  from  their  true 
values  agreed  well  with  the  estinrin  Si  i mrd  from  corresponding  accuracy 
studies  (Schwarz  ll)7f>).  11  should  * e noted,  however,  that  correlated  errors 

in  the  measurements  will  stronglv  influence  the  accuracy  of  the  results.  A 
second-order  Markov  sequence  wa-  u < d !>\  Schwarz  to  model  the  error  pro- 
cess along  the  profiles.  Depending  oi  the  bze  of  the  correlations  and  the 
variance  of  the  process,  the  mean  :-,|uar  rrors  would  tnoie  than  double  as 
compared  to  the  uncorrelated  ease,  .dm,1  url\,  a bias  in  the  data  would  impair 
the  accuracy  of  the  results  considerable. 

The  simulation  of  gravimetric  quantities  from  a point  mass  model  was 
considered  in  the  spectral  domain  and  conclusions  were  drawn  with  respect  to 
the  resulting  fields.  In  order  to  represent  adequately  regional  variations  of  the 
gravity  field  as  well  as  the  local  behaviour  of  the  second-order  gradients  the 
medium  and  the  high  frequency  part  of  the  spectrum  must  be  modeled  equally 
well.  This  can  only  he  achieved  by  vising  several  planes  of  generating  masses  at 
different  depths.  Schwarz  concluded  that  • many  cases  a model  with  two  planes 
may  already  l>e  sufficient. 
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2.3  Other  Stud  1 


I sing  the  results  i>f  the  ab..ve  studies,  Schwarz  computed  mean  gravity 
anomalies  from  g radiometer  and  altimeter  data  for  the  selected  test  area  which 
was  25u  in  latitude  times  : 0 m longitude.  The  results  of  these  extensive  com- 
putations as  well  as  the  eoi  r<  spending  eomputer  programs  w'ere  communicated 
to  A FGL. 

Schwa  tv  af  o did  tn  : . a eomldna  I 'on  of  satellite  derived  ha  rmonic  co- 
efficients and  terrestrial  men  cavity  anomalies  by  least-squares  collocation. 

The  results  were  communicated  informally  to  A FGL. 

3.  Satellite  to  Satellite  1* rack  ng  Rest  u h 

During  ih*  (ieriinl  of  the  e ntraet  an  evolutionary  research  development  has 
taken  plan  in  the  anahsi  of  saodlil*  I -.a  tell  i to  tracking  data  for  use  in  the  re- 
eovert  <d  m.  an  g i a v i ty  anomalie.  at  tlie  surface  of  the  earth.  This  work  is  de- 
scribed in  three  reports  winch  are  briefly  discussed  in  the  following. 

l lie  tirsi  report  in  this  an  a under  the  contract  was  that  of  Hummel,  Hajela, 
and  Hupp  , Hi  /(')).  fins  report  first  I'xamined  the  theory  where  one  satellite  is 
tracked  by  anotht  r such  that  a range  rate  lie  tween  a relay  satellite  (such  as  ATS-6) 
and  a close  satellite  (sueh  as  (tens- 3)  c tild  lie  determined.  This  range  rate  data 
can  be  used  to  determine  tlu  line  of  sight  acceleration  between  the  relay  and  the 
close  salt  llite  A number  of  different  geometries  were  considered  to  relate  this 
line  ot  sight  am  i ration  to  the  gradient  of  the  disturbing  potential.  For  one 
approximation  the  tol lowing  least  squa re  eollocation  solution  was  suggested; 

_ . •• RFS 

•^g(t^)  G(t  i)v.  , ( Ftr  jTr  j + ^ji ) ( li«c  cos  P ) i (19) 

where;  4g  is  the  predicted  anomaly; 

C((j  ) is  the  covariance  between  the  anomaly  being  predicted  and  the  ra- 
dial gravity  disturbance  at  the  close  satellite  observation  point  j; 

c.  i;  the  auto  covariance  matrix  between  the  radial  gravity  distur- 
bance components; 

I),.  is  the  noise  matrix 

UR  i.-  the  residual  line  of  - ighl  acceleration  with  respect  to  a refer- 
ence field  to  which  <Ag  will  lie  referred; 

R is  the  angle  between  the  line  of  sight  acceleration  at  the  close 
satellite  and  the  direction  of  the  radial  gravity  disturbance, 
bquution  ( ID)  i-  m approximation  in  the  sense  that  only  the  radial  component  of  the 
gravity  disturbance  is  being  considered. 

i lie  next  step  was  a series  of  numerical  simulation  studies  designed  to  test 
equation  (I'.b.  In  do  tins  orbits  were  generated  (using  Geodyne)  in  a reference 
field  and  in  a highei  degree  field  rlesig..  d to  reflect  reality.  Residuals  were 
formed  and  predictions  using  (1!))  were  carried  out  using  covariances  obtained 
from  ( OVAX  (Tseherning,  197(1).  Predictions  were  made  for  10°  and  5°  equal 
area  anomalies  with  different  data  point  intervals  and  different  values  for  the 
elements  of  the  I)  matrix,  in  this  work  it  was  estimated  that  with  the  geometry 
of  the  Deos-.l,  A IS  (.  satellite  situation,  1()'J  anomalies  could  be  recovered  to  an 
accuracy  (s.  d. ) of  about  + -I  mgals  and  a anomalies  to  an  accuracy  (s.d.  1 of  about 
± 1 1 mgals. 
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Additional  tests  indicated  that  the-  results  were  quite  sensitive  to  errors  in  the 
initial  orbital  elements  and  the  values  of  the  elements  in  the  I)  matrix. 

It  was  clear,  however,  that  the  suggested  procedure  did  work  in  a simula- 
tion mode  and  that  it  would  Lie  worthwhile  to  proceed  to  the  analysis  of  real  data. 

Such  an  analysis  is  described  in  the  report  of  llajela  (1977). 

In  the  research  carried  out  by  llajela  (ibid)  actual  Geos-3,  ATS-0  range 
rate  data  provided  by  NASA  was  analyzed.  The  data  supplied  was  in  the  ATS  It 
format  and  was  reformatted  for  entry  into  the  Geodyn  program.  Five  passes  of 
Geos-3  in  the  area  o : 15°  to  30°,  and  X:  27 D°  to  293°  wen?  selected  for  pro- 
cessing using  initial  state  vectors  supplied,  in  part,  by  NASA.  After  the  best 
set  of  initial  state  vectors  were  found  (w'ith  the  data  available)  the  range  rate 
residuals  were  fitted  to,  and  filtered  by  using  precise  continuous  cubic  splines 
to  fit  the  range-rate  data  in  the  least  squares  sense.  llajela  (ibid)  found  that  the 
optimum  spacing  of  the  nodes  where  the  adjacent  cubic  splines  met  was  even  (>0 
seconds  in  fitting  range-rate  observations  at  10  second  intervals.  After  the 
spline  fit  was  made,  the  residual  accelerations  ( it)  are  found  by  the  differenti- 
ation of  the  spline  function.  Using  this  process  a smoothed  set  of  accelerations 
were  obtained  for  parts  of  the  five  passes  of  Geos-3  across  the  area  of  interest. 

Hajela  (ibid)  used  this  data  to  rocmcr  light  f>°  equal  area  anomalies.  In 
doing  this  recovery  a number  of  different  variables  were  considered  including 
the  spacing  of  the  nodes  in  the  spline  fitting,  various  accuracy  estimates  for  use 
in  the  I)  matrix  of  equation  (1!)),  use  of  various  are  combinations,  and  consider- 
ation of  initial  epoch  vector  errors. 

The  recovered  anomalies  were  compared  to  ground  truth  data  available  from 
terrestrial  sources.  The  root  mean  square  discrepancy  between  the  recovered 
anomalies  and  the  ground  truth  values  was  H mgals  while  the  average  predicted 
standard  deviation  was  1 12  mgals.  It  was  felt  that  these  results  were  quite 
satisfactory  and  that  the  proposed  method  was  a workable  technique  for  gravity 
anomaly  recovery. 

After  the  above  study  it  became  clear  that  improved  modeling  could  be  done 
for  the  anomaly  recovery.  In  addition  a numlier  of  questions  remained  unanswered 
concerning  the  estimation  process.  The  main  problem  in  the  modelling  area  was 
in  relating  the  line  of  sight  accelerations  to  the  surface  gravity  anomalies.  In  the 
previous  reports  various  techniques  for  evaluating  the  covariances  were  described. 
However  the  application  was  not  as  direct  as  one  might  want.  To  improve  the  sit- 
uation Rummel  and  Rapp  (1977)  outlined  the  equations  needed  to  work  directly  with 
line  of  sight  accelerations  instead  of  radial  gravity  disturbances.  Hajela  (1978)  im- 
plemented these  equations  and  carried  out  new  tests  with  real  world  data  from  the 
ATS-6,  Geos-3  satellites. 

In  this  most  recent  study  Hajela  (ibid)  examined  a number  of  different  topics. 
One  case  considered  was  the  removal  of  a linear  trend  in  the  residual  accelerations 
caused  by  errors  in  the  initial  shite  vectors.  This  was  attempted  by  introducing  a 


systematic  part  (AX)  to  the  original  collocation  model.  Sample  solutions  were  made 


solving  for  a linear  trend.  However  Hajela  (ibid)  found  that  a considerable 
portion  of  the  signed  was  also  removed  from  the  residuals,  lhus  this  technique 
was  not  used  further. 

Hajela  also  looked  at  the  correlation  between  the  predicted  anomalies  from 
the  collocation  solution  using  the  data  from  the  previous  solution.  He  first  made 
ti  sts  for  the  recovery  of  the  5°  equal  area  anomalies.  He  found  for  the  eight 
adjacent  anomalies  the  largest  correlation  coefficient  was  -0.09  indicating  neg- 
ligible correlation.  A second  test  was  made  for  the  recovery  of  2.°5  anomalies. 
The  largest  correlation  found  for  a 30  second  data  interval  was  0.1  which  increased 
to  0.3  when  a 4 minute  interval  was  used.  The  biggest  change  (from  the  5r'  solu- 
tion) was  an  increase  in  the  standard  deviation  from  ±10  mgals  for  the  5°  anom- 
alies to  <- 15  mgals  for  the  2.°  5 anomalies. 

I ! c main  purpose  of  the  Hajela  (1978)  report  was  the  implementation  of  the 
line  of  sight  acceleration  method.  Hajela  gives  the  specific  equations  for  the 
rigorous  implementation  of  this  method.  He  then  uses  this  method  with  real  data 
based  on  improved  satellite  orbits.  Data  from  a number  of  arcs  was  provided  to 
us  by  Jim  Marsh  from  NASA.  Hajela  used  the  new  technique  for  5°  anomaly  re- 
covci  > „ The  recovered  anomalies  ineight  blocks  were  compared  to  corresponding 
terrestrial  anomalies  and  anomalies  recently  derived  from  Geos-3  altimeter  data. 
The  root  mean  square  difference  between  the  SST derived  anomaly  and  the  altimctei 
derived  anomaly  was  ±7  mgals  with  an  estimated  standard  deviation  of  each  SST 
anomaly  being  approximately  ± G mgals. 

These  new  tests  indicate  that  the1  proposed  method  is  a valid  technique  for 
anomaly  recovery  from  SST  data.  We  now  need  more  such  data  for  further  testing 
of  the  method. 


4 . c onve  r go  nee  Problems 

Despite  the  general  feeling  that  the  convergence  of  the  spherical  harmonic 
expansions  to  finite  degree,  on  the  surface  of  the  earth,  is  not  a problem,  we  have 
carried  out  several  studies  to  improve  our  understanding  of  the  problem. 

One  such  study,  Rapp  (1977b),  described  and  tested  a method  for  computing 
gravity  anomalies  at  the  surface  of  the  earth  from  spherical  harmonic  expansion 
in  such  a way  ns  to  avoid  convergence  questions.  The  method  is  simple, being 
basically  a two  step  process.  The  first  step  takes  the  spherical  harmonic  ex- 
pansion of  the  disturbing  potential  and  then  evaluates  it  at  points,  or  in  compoart- 
ments,  on  a sphere  that  surrounds  the  mass  of  the  earth.  (The  mass  of  the  atmos- 
phere was  not  considered  in  these  computations.)  The  anomalies  are  then  down- 
ward continued  to  the  surface  of  the  earth  using  least  squares  collocation  techniques. 
To  test  this  idea,  f>°  anomalies  were  used  with  their  value  computed  at  the  surface 
of  the  earth  in  two  ways.  First  the  anomalies  were  computed  on  the  bounding 
sphere.  A downward  continuation  correction  was  applied  to  obtain  a surface  anomaly^ 
that  was  corrqwi red  to  the  terrestrial  data  where  a mean  square  difference  of  91  mgal 
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was  found.  Then  the  potential  coelT.eienl  anomalies  were  eomputed  dirccth  ,1 
the  surface  where  a mean  square  difference  of  10!)  mgal  was  found.  In  this 
case  the  theoretically  more  correct  procedure  gave  the  heller  result. 

The  concept  of  the  enclosing  sphere  computation  also  provides  a convenient 
formulation  of  the  problem  of  computm-  potential  coefficients  from  surfaci  cav- 
ity anomalies.  This  procedure  can  ■ ■ r,  presented  in  the  following  form: 

('  <'»  + £<'  * A<\,  (20) 


where:  C is  the  rigorous  potential  e ■ fficient; 

CA  is  the  approximate'  potent : ' coefficient  derived  from  the  usual  sum 

mation  formula  applied  to  uncorrccted  anomaly  data; 

AC|  is  a correction  dependenf  on  the  radius  of  the  hounding  sphere,  a 
mean  earth  radius,  th < equal  trial  radius  and  th<  (',  values; 

A Ci i is  a correction  term  computed  by  applying  the  usual  summation 
formula  to  the  upward  continuation  correction  term. 

Studies  were  made  of  the  magnitude  of  Ac,  and  expressed  as  a percentage 

of  the  expected  magnitudes  of  tin  haul  coefficients.  Up  to  dcgrei  12  the  total 
correction  was  less  than  2',  . At  d<  groo  lo  the  total  correction  had  reached  7 . 
These  figures  indicated  the  corrections  t>  rms  are  small  and  negligible  with  the 
current  data  accuracy.  Additional  invi  stigalions  to  highei*  degrees  are  needed 
using  1°  x 1°  anomalies. 

A completely  different  view  of  this  problem  was  taken  by  Sjoberg  (1977a). 

In  the  first  part  of  his  study  Sjoberg  const  meted  a simple  example  to  show  that 
the  spherical  harmonic  expansion  of  th,  > arth's  gravitational  potential  is  diver- 
gent at  the  surface  of  the  earth.  Th,  < r r caused  |,y  evaluating  a divergent 
series  but  including  only  a finite  number  of  terms  was  considered  in  terms  of  an 
analytic  continuation  error  and  a truncation  error.  The  percentage  error  in  the 
disturbing  potential  at  degree  (it)  could  reach  207,.  when  an  improper  downward 
continuation  procedure  was  applied  to  a sphere  with  a point  mass  located  at  a 
height  of  20  km  above  the  sphere. 

Sjoberg  (ibid)  extended  his  special  case  to  that  of  the  real  earth  containing 
topography.  He  computed  errors  in  the  potential  caused  by  neglecting  topography 
and  compared  his  results  to  previous  results  found  by  (,'ook  and  Uevallois.  Specific 
equations  for  computing  corrections  term  due  to  the  topography  were  derived  and 
applied  for  a spherical  harmonic  expansion  of  the  disturbing  potential  and  the  grav- 
ity anomaly.  These  equations  were  applied  to  obtain  global  error  estimates  and  to 
obtain  specific  error  estimates  using  5°  x .7°  mean  elevations.  Krrors  on  the  order 
of  ±1  meter  were  found  for  geoid  undulations  and  up  to  70  mgals  for  the  anomalies. 
This  large  anomaly  error  is  unrealistic  so  that  there  is  some  concern  on  the  cor- 
rectness of  some  of  the  equations  or  data  used  by  Sjoberg  (ibid).  It  may  be  that 
smaller  blocks  (say  1°  x lc  ) are  needed  to  adequately  represent  the  true  situation. 
Thus  additional  work  is  called  for  in  this  area. 


• >.  ( Hhcr  A i t :is  of  Resca  rch 

\s  pari  if  on  r continuing  into  rest  in  estimating  accuracies  of  gravimetric 
quantities  SjdinTg  (IU77b)  ( (inside iv<l  (In  estimation  of  the  accuracies  of  the 
computation  <>1  lli«  deflection  "f  the  vertical.  In  this  estimation  Sjoberg  (ibirl) 
considered  that  (lie  deflections  would  lie  computed  in  three  components:  1 ) a 
e.  nlriliulion  from  a set  of  potent  ml  eoeffii  ients;  2 ) a contribution  from  l°x  1° 
uion.alies  in  an  area  u rmiindm;:  die  eom()utation  point  ; and  .‘5  ) a contribution 
from  detailed  data  in  the  immediate  vicinity  of  the  comjiutation  point.  Since  po- 

• ential  e .1  ffiii.nl.  ire  Known  t • a relatively  low  degree  it  was  also  necessary 
t"  i on  ider  a trim  at  ion  error.  Also,  the  use  of  1°  x 1°  anomalies  implies  that 
i * rlain  higher  In  e u ncy  information  will  lie  lost  in  the  solution  so  that  this 
effect  must  be  Considered. 

for  each  ot  111.  •se  cffci  Is  nr  contributions,  Sjoberg  (ibid)  derived  approx- 
imaie  error  equations.  lie  applied  these  equations  with  the  OKM  7 potential  co- 
eflicients  and  existing  f x 1°  mean  anomaly  information.  Specific  accuracy 
estimate:  due  ;.>  the  \arious  error  sources  were  made  for  four  points  having 
various  accuracies  of  1 x 1°  anomalies  surrounding  them.  Considering  all 
error  sources  w itli  tin  local  data  surrounding  the  point  taken  out  to  1°  from  the 
point,  Sjolx  rg  (ibid)  estimated  the  error  in  the  total  deflection  to  range  from 
± 1.  > ’ to  +.2.  .i"  with  a cap  to  20f  of  surrounding  1°  x 1°  data. 

In  another  a r<  a Kearslev  ( 11)77)  examined  the  prediction  of  mean  anomalies 
a!  sea  using  collocation  techniques  applied  to  selected  geophysical  phenomena.  In 
this  paper  the  author  reviewed  some  of  the  factors  that  influence  anomalies  at  sea 
including  bottom  topography  (or  ocean  depth),  and  the  age  of  the  crust  in  the  area 
of  interest. 

In  considering  these  dependencies,  Kearslev  (ibid)  developed  auto  covariance 
tune  lions  for  crustal  age,  the  re-:  covariance  function  for  anomalies  and  crustal 
age,  auto  ovarian'  ■ functions  for  depth  and  a cross  covariance  function  for  anom- 
alies and  depth,  all  using  data  in  1“  x V blocks.  This  information  was  used  to 
predict  a number  of  1 x V'  anomalies  which  were  then  compared  to  their  known 
values.  The  root  mean  square  difference  found  was  +.11.5  mgals  which  compares 
favorably  with  the  predicted  standard  deviations  of  about  +11  mgals.  Kearslev 
(il.id)  suggests  that  all  the  data  and  covariances  could  be  combined  in  one  general 
prediction  process. 

in  (In  |iast  few  years  great  progress  has  been  made  in  improving  our  knowl- 
edge of  the  earth’s  gravity  field,  potential  coefficient  determinations  from  sat- 
ellite data  alone  have  Ixu  n extended  to  degree  20  with  many  additional  resonant 
type  terms.  In  a number  of  new  areas  terrestrial  data  has  become  available. 

And  of  most  significant  importance  the  Geos-.')  altimeter  data  has  enabled  the 
determination  of  a large  number  of  1°  x l°  mean  anomalies.  Rapp  (1978)  suggests 
and  implements  a procedure  to  combine  all  this  data  in  a rigorous  least  squares 
adjustiiK  nl  process.  1’h.  technique  was  applied  to  the  GEM  9 potential  coefficient 
set  taken  to  degree  1 2 field  on  an  a linos'  n set  of  1°  x 1°  anomalies.  Speci- 
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fically  this  anomaly  field  contained  50050  values,  2s|7<i  of  which  w.  n derived 
from  Geos-3  altimeter  data.  The  remaining  values  needed  ti>  form  a g! . t,:d  s<  i 
were  taken  to  be  zero  with  a standard  deviation  of  > :to  mgals.  The  result  of  the 
adjustment  was  an  adjusted  set  of  |M>tential  coefficients  to  degree  12  and  an  ad- 
justed global  set  of  1 x I anomalies.  These  adjusted  anomalies  were  deve|o|>ed 
into  potential  coefficients  to  degree  00.  These  coefficients  must  agree  exartlx 
with  the  coefficients  found  in  the  actual  adjustment  process.  Higher  degree  so 
lutions  are  possible,  both  in  the  adjustment  process,  and  in  the  determination  of 
the  potential  coefficients  from  the  adjusted  anomalies.  The  main  advantage  of  the 
technique  used  in  liapp  (ibid)  is  that  it  provides  a rigorous  combination  of  the  data 
to  obtain  a consistent  set  of  potential  coefficients  and  1°  x 1°  anomalies. 


(i. 


References 


<>.  1 References  i‘ reduced  Under  the  ( 'ontract 

Unjoin,  |).  p, , Recovery  of  5”  Mean  Gravity  Anomalies  in  Local  Areas  From 
ATS-6.  GKON-3  Satellite  to  Satellite  Range-Rate  Observations,  Report 
No.  259,  Department  of  Geodetic  Science,  The  Ohio  State  University, 
AFGL  l'R-77-0272,  Scientific  Report  No.  11,  September  1977. 

Ilajela,  1).  P. , Improved  Procedures  for  the  Recovery  of  5°  Mean  Gravity 
Anomalies  From  ATS-6/GKOS-3  Satellite  to  Satellite  Itange-Rate 
Observations  Using  Least  Squares  Collocation,  Report  No.  276, 
Department  of  Geodetic  Science,  The  Ohio  State  University,  AFGL-TR- 
79-0260,  Scientific  Report  No.  20,  September  1978. 

Jekcli,  Christopher,  An  Investigation  of  Two  Models  for  the  Degree  Variances 
of  Global  Covariance  Functions,  (Thesis),  Report  No.  275,  Department 
of  Geodetic  Science,  The  Ohio  State  University,  AFGL-TR-78-0235,  Sci  - 
entific Report  No.  IS,  September  1978. 


Kearsley,  William,  The  Kstimation  of  Mean  Gravity  Anomalies  at  Sea  from 

Other  Geophysical  Phenomena,  Report  No.  270,  Department  of  Geodetic 
Science,  The  Ohio  State  University,  A FGL-TR-78-0069,  Scientific  Report 
No.  14,  December  1977. 

Moritz,  Helmut,  Integral  Formulas  and  Collocation,  Report  No.  234,  Department 
of  Geodetic  Science,  The  Ohio  State  University,  AFCRL-TR-75-0628, 
Scientific  Report  No.  1,  December  1975. 

Moritz,  Helmut,  Covariance  Functions  in  Least -Squares  Collocation,  Report  No. 
240,  Department  of  Geodetic  Science,  The  Ohio  State  University,  AFGL- 
TR-7G-0165,  Scientific  Report  No.  3,  June  1976a. 

Moritz,  Helmut,  Least-Squares  Collocation  as  a Gravitational  inverse  Problem, 
Report  No.  249,  Department  of  Geodetic  Science,  The  Ohio  State  University, 
A FG l j-T R- 7 6-0 278,  Scientific  Report  No.  6,  November  1976b. 

Moritz,  Helmut,  On  the  Computation  of  a Global  Covariance  Model,  Report  No. 
255,  Department  of  Geodetic  Science,  The  Ohio  State  University,  AFGL- 
TR-77-0 163,  Scientific  Report  No.  10,  June  1977. 

Moritz,  Helmut,  Statistical  Foundations  of  Collocation,  Report  No.  272,  Depart- 
ment of  Geodetic  Science,  The  Ohio  State  University,  AFGL-TR-78-0182, 
Scientific  Report  No.  16,  June  1978a. 

Moritz,  Helmut,  The  Operational  Approach  to  Physical  Geodesy,  Report  No.  277, 
Department  of  Geodetic  Science,  The  Ohio  State  University,  A FGL-TR-78- 
0281,  Scientific  Report  No.  21,  October  1978b. 

-26- 


Rapp,  R.  H. , The  Use  of  Gravit\  Anomalies  on  a Mounding  Sphe re  to  Improve 
Potential  Coefficient  Determinations,  Report  No.  251,  Department  of 
Geodetic  Science,  The  ( ihio  Slab  University,  Al'Gl,  TIG77  OH(i, 

Scientific  Re|jort  No.  t) , June  10,71,. 

Rapp,  It.  H.  , A Global  1 v 1 Anomalv  Field  Combining  Satellite,  (IKOS-l! 
Altimeter  and  Terrestrial  AnomaK  Data,  Iteport  No.  278,  Di'partment 
of  Geodetic  Science,  rhe  ( i'h,,  sbio  University,  Al'Gl,  TR-7H-0282, 
Scientific  Iteport  No.  72,  Septemlier  lays. 

Rummel,  Reiner,  A Model  Compa  ris or  , i ■ast  -Sepia  res  ( ollocat  ion.  Report 
No.  238,  Department  of  Geodeli  Sou  nee,  The  Ohio  State  UniversiU, 

A FG  1,-1  11—7 G — 0 19 8 , Scientifu  Rep,  ,rl  No.  5,  August  1070. 

Rummel,  Reiner,  D.  P.  Hajela,  and  R.  |i.  Rap|,f  Recove r\  of  Mean  Gravity 

Anomalies  from  Satellite-Satellite  Range-Rate  Data  Using  feast  Squares 
Collocation,  Report  No.  2f  , f,  i arlmenl  of  Geodetic  Science,  The  Ohio 
State  University,  AFGf-IU  (o  0201,  Scientific  Report  No.  7,  September 
1976. 

Schwarz,  Klaus -Peter,  Geodetic  Accuracies  obtainable  from  Measurements  of 
First  and  Second  Order  Gravitational  Gradients,  Report  No.  212,  Depart- 
ment of  Geodetic  Science,  The  Ohio  State'  University,  A FGf-TR-76-Ols!), 
Scientific  Report  No.  1,  August  1976. 

Schwarz,  Klaus-Peter,  Simulation  stuff  of  Airborne*  Gradiomelry,  Report  No. 

253,  Department  of  Geodetic  Science.  The  Ohio  Slate  University,  Al'Gl  I'll 
77-0129,  Scientific  Report  No.  s,  ypiv  1977. 

Sjoberg,  Iars,  On  the  Errors  of  Spherical  Harmonic  Developments  of  GraviU  at 
the  Surface  of  the  Earth,  Re  port  No.  257,  Department  of  Geodetic  Science, 
The  Ohio  Shite  Universitv,  Al'Gl.  IR-77-0229,  Scientific  Report  No.  12, 
August  1977a. 

Sjoberg,  lars,  The  Accuracy  of  Gra\  air  I ric  Deflections  of  the  Vertical  as 
Derived  from  the  Gl-'M  7 Potential  Coefficients  and  Terrestrial  Gra\it\ 

Data,  Report  No.  265,  Department  of  Geodetic  Science,  The  Ohio  Stub 
University,  A FGf-T  R-77-02H7,  Scientific  Report  No.  13,  Novcmlnr  I977l>. 

Sjoberg,  Lars,  Potential  Coefficient  Determinations  from  10°  Terrestrial 
Gravity  Data  by  Means  of  Collocation,  Report  No.  27-1,  Department  of 
Geodetic  Science,  The  Ohio  State  I diversity,  A FGI.-TR-78-02I1 , Scientific 
Report  No.  19,  September  1978a. 

Sjoberg,  Uirs,  A Comparison  of  Rjerham mar’s  Method  and  Collocation  in 

Physical  Geodesy,  Report  No.  273,  I )epa  rtment  of  Geodetic  Science,  The 
Ohio  State  University,  Al'Gl,  TR-78-0203,  Scientific  Ri']»ort  No.  17, 
Septcmlx'r  1978b. 


Siinkcl,  Hans,  Approximation  of  Covariance  Functions  by  Non-Positive  Functions, 
Report  No,  tl 7 l , Department  of  Geodetic  Science,  The  Ohio  State  University, 
A Ft! I .-TR-78-0  177,  Seientifie  Report  No.  l.r>,  May  1978. 

I'seherning,  C.C.,  C'ovarianee  Kxpressions  for  Second  and  Lower  Order  Deriva- 
tives of  the  Anomalous  Potential,  Repoi’t  No.  225,  Department  of  Geodetic 
Science,  The  Ohio  State  University,  AFC.B-TR-76-0052,  Scientific  Report 
No.  2,  January  1970. 

<>.  2 ( )ther  References 


Ruckus,  G. , inference  from  Inadequate  and  Inaccurate  Data,  land  II,  Proc. 

Natl.  Acad.  Sci.,  (15,  1-7  and  281-287,  1970. 

Krarup,  I.,  A Contribution  to  the  Mathematical  Foundation  of  Physical 
Geodesy,  Publ.  41,  Danish  Geodetic  Institute,  Copenhagen,  1969. 

I,:iu  rit/.en,  S.l,.,  The  Probabilistic  Background  of  Some  Statistical  Methods  in 
Physical  Geodesy,  Publ.  No.  48,  Danish  Geodetic  Institue,  Copenhagen, 
1972. 

Meiss],  P. , Probabilistic  l-.rror  Analysis  of  Airborne  Gravimetry,  Report  No. 
128,  Department  of  Geodetic  Science,  The  Ohio  State  University,  Columbus 
1970. 


Moritz,  ||.,  Kinematical  Geodesy,  Report  No.  92,  Department  of  Geodetic  Science, 
The  Ohio  State  University,  Columbus,  1967. 

Moritz,  ||.,  Advanced  Beast-Squares  Methods,  Report  No.  175,  Department  of 
Geodelie  Science,  The  Ohio  State  University,  Columbus,  1972. 

'I'Til.  , II.,  and  K.  P.  Schwarz,  On  the  Computation  of  Spherical  Harmonics  from 
Satellite  observations,  Rulletinn  di  Geodcsia  e Scienze  Affini,  32,  3,  185- 
200,  1972. 

Moritz,  ||.,  Precise  Gravimetric  Geodesy,  Report  No.  219,  Department  of 
Geodetic  Science,  The  Ohio  State  University,  Columbus,  1974. 

Rapp,  R.  II.,  Potential  Coefficient  Determination  from  5°  Terrestrial  Gravity 
Data,  Report  No.  251,  Department  of  Geodetic  Science,  The  Ohio  State 
University,  Columbus,  1977a.  Jk. 

Sza bo.  It.  and  |).  Anthony,  Results  of  A FCRB’S  Kxperimental  Aerial  Gravity 
Measurements,  Bulletin  < .codesiquc,  No.  100,  1971. 

I scheming,  C.c,  and  R . II.  Rapp,  Closed  ('ova ria nee  Kxpressions  for  Gravity 
Anomalies,  Geoid  Undulations  and  Deflections  of  the  Vertical  Implied  by 
Anomaly  Degree  Variance  Models,  Report  No.  208,  Department  of  Geodetic 
Science,  The  Ohio  Stale  University,  Columbus,  1974. 

-28-' 


7.  List  of  Scientific  IVrsonni‘1 


Supervisors  and  Principal  Investigators 

Urho  A.  Uotila 
Richard  H.  Rapp 


Research  Associates  and  Consultants 

D.  P.  Hajela 
William  Hears  lev 
Helmut  Moritz 
Reiner  Rummel 
Klaus-Peter  Schwarz 
lars  Sjoberg 
Hans  Siinkcl 
C.C.  Tscherning 


Graduate  Research  Associates 

Joshua  Greenfield 
Christopher  Jekeli 
D.  Jeyanandan 
Lenny  A.  Krieg 
Joseph  C.  Loon 
K.  Katsamhalos 


